Fluctuation relations with intermittent non-Gaussian variables 
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Non-equilibrium stationary fluctuations may exhibit a special symmetry called fluctuation rela- 
tions (FR). Here, we show that this property is always satisfied by the subtraction of two random 
and independent variables related by a thermodynamic-like change of measure. Taking one of them 
as a modulated Poisson process, it is demonstrated that intermittence and FR are compatibles 
properties that may coexist naturally. Strong non-Gaussian features characterize the probability 
distribution and its generating function. Their associated large deviation functions (LDF) develop 
a "kink" at the origin and a plateau regime respectively. Application of this model in different 
stationary nonequilibrium situations is discussed. 
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I. INTRODUCTION 

Intermittency is a phenomenon that appears in a wide 
class of situations such as chaotic dynamics [l| , atomic [1| 
and nanoscopic Q fluorescent systems, single-molecule 
reaction dynamics [3], biological self-organized models 
or fluid turbulence just to name a few. It consists 
in a random switching amongst qualitatively different 
system dynamical regimes. Its stochastic behavior may 
be ergodic or not. Complexity, non-Gaussian statistics, 
and nonequilibrium dynamics are closely related with its 
development. In this last context, (Gallavoti-Cohen) FR 
[tII addresses the distribution of entropy production in 
far-from-cquilibrium steady states. It relates the proba- 
bility of observing a certain entropy production rate to 
the probability of observing the corresponding entropy 
consumption rate 0412 • Theoretical and experimental 
results confirm its validity in diverse nonequilibrium sys- 
tems [T3l - [23 |. Interestingly, the symmetry imposed by 
FR has also been found valid for non-thermodynamic 
variables [23} . 

As FR and intermittence are intrinsically related with 
nonequilibrium dynamics, it is natural to ask about the 
possible coexistence of both properties. The main goal of 
this contribution is to give a positive answer to this issue. 
We demonstrate that there may exist non-equilibrium 
steady states whose (non-Gaussian) fluctuations are in- 
termittent and satisfy the FR. An explicit construction of 
a stochastic variable with the required symmetry defines 
the basis of the present analysis. It is shown that FR can 
always be satisfied by the subtraction of two independent 
stochastic variables related by an exponential change of 
measure. Hence, intermittence is introduced by choos- 
ing one of them as a random modulated Poisson process. 
Special attention is paid to the asymptotic regime, where 
FR can be analyzed through the large deviation functions 
(LDF) 0-11] of the probability distribution and its associ- 
ated generating function. In the parameter regime where 
intermittence arises, they develop a kink around the ori- 
gin and a plateau regime respectively. These strong non- 
Gaussian features are similar to that found in different 



observables such as the entropy production rate of a col- 
loidal particle [22i] and the (time-average) velocity of a 
polar granular rod [l^. While the specific analysis of 
these systems is beyond the scope of this contribution, 
the present results suggest that intermittence under the 
constraint of FR symmetries may be a central ingredi- 
ent when studying nonequilibirum states characterized 
by non-Gaussian fluctuations. 



II. FLUCTUATION RELATIONS WITH 
INDEPENDENT STOCHASTIC VARIABLES 

Our results rely in the following analysis. Let an arbi- 
trary stochastic variable Xst, whose probability distribu- 
tion P{x) satisfies a FR defined as 



P{~x) = P(a;)exp[-Ca;], 



(1) 



where is a given real positive constant. This symme- 
try implies that negative values are exponentially less 
probable than positive values. In terms of the generating 
function Z{X) = (e"^^"') = dxP{x)e-^'' , Z{0) = 1, 
the FR reads 



Z(A) = Z(-A + C). 



(2) 



In order to satisfy the FR, we write Xst as the subtraction 
of two statistically independent variables xf^ , 



Xst 



(3) 



Then, Z(A) = Z+(A)Z_(-A), where Z±(A) are the gener- 
ating functions of x^^., and whose distributions are P±{x). 
The condition ([2]) is satisfied by demanding the relation 



Z_(A) 
which in turn implies 



z+jx + O 
z+iO ' 



P.ix)=P+{x) 



(4) 



(5) 
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where (e ''^)^ = Z+{C). This is one of the central results 
of this paper. It says us that the FR symmetry ([T|) is 
always satisfied by ^ whenever the change of measure 
([5]) is imposed. By reading x as an energy and C as an 
inverse temperature, its thermodynamic like-structure is 
self-evident. Similar transformations were introduced in 
Refs. [2i - l26| . The structure of Eq. ^ can also be read 
as a variation of the change of measure used in umbrella 
sampling [27j . 

Note that in Eq. ([5]) no conditions are required on 
P+{x), while P-{x) is determined from it. Alternatively, 
one can choose P_ (x) as the independent distribution. 
In fact, Eq. ([2]) can also be satisfied by taking 



Z-{x-0 
Z-{-0 ' 



delivering the relation 



P+{x)=P-{x) 



(6) 



(7) 



where ^e'*'''^)_ = C). These last two expressions 

follow straightforwardly by inverting Eqs. ^ and ([5]). 
Nevertheless, what is important to realize is that, in gen- 
eral, both kind of solutions lead to different family of 
solutions if one assume the same property for each "in- 
dependent" distribution. 

The relations (O and (O do not involves time. There- 
fore, the proposal Eq. ([3]) guarantees the fulfillment of 
the symmetry ([1} at any time. 



Time-average variables 

When studying a given nonequilibrium steady state, 
the stochastic variable of interest may be defined through 
a time average, 
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t{v) 



*^^dt',^,(t')^^£^(i±4_^, (8) 

t{v) 



where [d/ dt)xst{t) = Vst{t). This definition makes sense 
when the average of the integral contribution grows lin- 
early with time. For example, Xst{t) may represents en- 
tropy production and Vstit) its rate [22], or respecti vely 
the position and velocity of a self-propelled particle [23 . 
{v) is the stationary mean value of Vst(t). 

Usually the regime of interest is the stationary one, 
t = oo. Hence, we write fj,r = limt^cxj MT(i)- For this 
stationary variable, instead of Eq. ((T)) , we write the FR 



P{-I^t) = P{^J.T)cxp[-a^J.rT]■ 



(9) 



This condition can also be satisfied in the present ap- 
proach because the statistical properties of follows 
from those of the stationary increments of Xst{t), that 
is x^ = \imt^oo[xst{t + t) - Xst{t)]. By writing x^ 



as the subtraction of two independent stochastic vari- 
ables related by the change of measures ([5]) [or ([7])], it 
follows P{—x^) = exp[—(x^]P{x'^). The connection 
between x'^ and fir implies the probabilities relation 
P{fir)dfir = Pix^)dx^ [28], which lead to FR Eq. ^ 
with a = C (u) . Consequently, the proposal based on in- 
dependent stochastic contributions, as in the previous 
case [Eq. ([Ij], allows to fulfill the symmetry © for any 
value of r. 



Large deviation functions 

In general, a (rate) variable [Eq. ([5])] that characterizes 
a given non-equilibrium state only satisfies the relation 
© in a long time regime, r — >■ oo. In fact, the Gallavoti- 
Cohen FR is an asymptotic relation in time 0-11] • In 
this regime, when the probability distribution adopt the 
asymptotic structure (limT-_j.oo) P{^-t) ~ exp[— T(^(/iT-)], 
the statistics can be analyzed through a large deviation 
theory [2^. Consistently, the generating function Z{X) = 
^^-Xth^'^ scales in the same way, Z{X) « exp[— r0(A)]. 
Both ifiiidr) and 6(A) define the LDF's of the problem. 
They completely characterize the asymptotic regime. In 
terms of (fifJ-r), the FR reads 



1 



■log 



-ip{+flr) + Vi-fJ-r) ^ O-I-I'T- (10) 



Through a sadle-point approximation, both LDF can be 
related by a Legendre-Fenchel transformation 24] 

93(/ir) = max[0(A) — A/ir], = min[93(/ir) + A/ir]. 

(11) 

These relations and Eq. ((TU| lead to the equivalent for- 
mulation of the FR symmetry [9|] 

e(A) = e(-A-l-a). (12) 

Conditions (|T0| and (IT2]) are well known expressions of 
Gallavoti-Cohen FR symmetry 0-[l]. As the independent 
variables splitting [Eq. (O] allows to fulfill the FR (O at 
any time, trivially it can also be utilized in the long time 
regime. After explicitly writing fir = /i^ — /-i^, Eq. (|4]) 
and the asymptotic structure of Z{\) allows us to write 

e(A) = e+(A) + e+(-A + a)-e+(a), (13) 

where 0+(A) is the LDF corresponding to the generat- 
ing function of /x^. By knowing (an arbitrary) 0+(A) the 
previous expression defines 6(A), which in turn through 
Eq. (fTTj) provides the LDF (fljir). By construction, the 
fulfillment of conditions pH)) and ([T^ is guaranteed. Al- 
ternatively, from Eq. (jS]) we can also write 



e(A) = e_(-A) + e_(A - a) - e_(-a). 



(14) 



where now 6- (A) is the LDF associated to the gener- 
ating function of . The relations ([13]) and ([Ti]) define 
the second main result of this paper. They allow to char- 
acterizing the long time regime in terms of asymptotic 
properties of the independent stochastic contributions. 



3 



III. INTERMITTENT VARIABLES 

Our results allows to build up a variable that satisfy 
the FR Eq. ([9|) after knowing the statistical properties 
of an arbitrary one. The regime of interest is the sta- 
tionary one, where the FR symmetry is characterized by 
Eq. ([T3|) or (fT4|) . For example one can assume Gaus- 
sian or Poissonian statistics. In this last case, if one take 
0+(A) = 7(1 — e^'^), which correspond to the LDF of an 
unidirectional Poisson process with rate 7, from Eq. (I13|) 
we get e(A) = 7(1 - e"^) + 7e""(l - e-^). We note that 
this expression corresponds to the LDF for the entropy 
production of an asymmetric random walk [22|. Hence, 
while the previous analysis seem to be rather abstract, 
there exist non-trivial dynamics where they apply. 

The random walk model is able to fit some non- 
Gaussian properties found in Ref. I2a. Here, motivated 
by the experimental results of Ref. |23| . we introduce a 
similar generalized model able to develop intermittence. 
On the basis of previous analysis, the velocity [see Eq. 
([H])] is defined as Vstit) — — v~i{t), where 



(15) 



The constant xq introduces the right units of Vstit), 6{t) 
is the Dirac delta function, and tf are successive random 
times. Consistently, Xst{t) — x'^ti^) ~ 



(16) 



where nff{t) is the number of (idelta) events in the 
interval (0,t). Hence, they are positive (discrete) ran- 
dom variables. Their generating functions are denoted 
as Z±(s,0 = (e— ^ E™=o9;^We-^", where 
{9m(0}m=o the respective counting probabilities. In- 
dependently of their statistics, by choosing n^j(t) as the 
"free" variable, this model satisfy the FR symmetry (|9]) 
after demanding [see Eq. ([S])] the change of measures 



^e-son^t(t)) ' 



(17) 



where (e"''''"^*'*-') = Z^{so, t), and sq is an arbitrary pos- 
itive dimensionless constant. Hence, njf{t) can be read 
as the "sq— ensemble" associated to n'^th') ^6]- After a 



change of variables based on Eqs. ([5]) and p^ . the con- 
stant a reads 

a = soSI = so(/+ - /-), (18) 

where = —\imt^oot~^{d/ds)Zi^{s,t)\s=o, or equiv- 
alently liijit-foo {nff{t)) ~ I^t, which in turn implies 
{v} ^xo{I+ 

In order to close the model, it is necessary to specify 
the statistical properties of (t) . Its generating function 
is written as 
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FIG. 1: (color online) (a) Probability distribution P{i-Lt) that 
follows from Eq. (|20p (see text) for different values of r. (b) 
Linear dependence of (1/r) In [P(-|-/iT-)/^'(— Mr)] correspond- 
ing to the distributions shown in (a). The value of the slope 
a was obtained from Eq. pSI) . The parameters [Eq. (|20p ] are 
so = 0.2, 7 = WOs-\ Ta = 0.54s-\ and Lj = 4s-\ 



where the evolution of each contribution read 



dt 

dZj{s,t) 
Jt 



= -esZ+{s,t) ~VAZ\{s,t) + ViZ+{s,t), 

= +TAZ+{s,t)-TiZ+{s,t), (20) 



with 9s = 7(1 



'). These dynamics allow us to read 



Z+(s,t) =Z+(s,t) + Z+(s,t), 



(19) 



( {t) as a modulated Poissonian (counting) process [2£ 
whose rate at random times adopts the values 7a = 7 
(Active regime) and 7/ = (Inactive regime). The 
switching between both states is governed by a classical 
master equation with transition rates F^ and F/. The 
asymmetric random walk model 2^] is recovered from 
Eq. (gni) by taking F^// = 0. 

The generating function Z~{s, t) can be obtained from 
Eq. Q after knowing Z+(s,i). Hence, the (lattice) dis- 
tribution [2^ of the process Xst{t) can be obtained by 
finding Z^{s,t) from Eq. (1201) and a posterior numerical 
Fourier inversion of Zn{s,t) — Z^ {s,t)Z^ {~s,t) in the 
s — —ik variable. By using the Markovian property of 
Eq. ((20)) . the probability of x"^ = liuit^oo Xst{t) follows 
from that of Xstit) by taking stationary initial conditions 
for the rate fluctuations. The distribution P{^t) of Mr 
follows from the change of variables defined by Eq. ([5]) . 

In Fig. 1(a) we plot the distribution Pinr) for dif- 
ferent values of t. For all times the distributions satisfy 
the symmetry Eq. (jH]) [Fig. 1(b)]. As expected, we con- 
firmed that this property is valid for any value of the 
rate parameters that define the evolution (|20| . On the 
other hand, the short time behavior of PifJ-r) strongly 
depends on the chosen parameters values. Nevertheless, 
for increasing r the distributions becomes similar and, 
consistently, develop an increasing peak around fXr — 
In this regime [r > 0.3 in Fig. 1(a)] a large deviation the- 
ory applies. Hence the problem can be analyzed through 
the corresponding LDF. 
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Large deviations functions 

The long time behavior of Z^{s,t) completely define 
the asymptotic statistical properties of the counting pro- 
cess n^ti^)- Before characterizing the LDF functions from 
it, we notice that in the long time limit some results 
can be established for Z~{s,t). As demonstrated in Ref. 
[26l |. n'^t{t) can be mapped with a renewal process char- 
acterized by a shift closure property, which implies that 
asymptotically n~^(t) also becomes a renewal process 
with a renormalized waiting time distribution. Hence, 
asymptotically Z~{s,t) has the same structure and dy- 
namics than Z^{s,t) but with renormalized rates. After 
some hard calculations steps, which are not relevant for 
the following analysis, we obtained 7 -> jsg ~ 7e^^", 
and more complex expressions for the hopping rates 
^A/i ^ TA/iiso). Furthermore /+ = iTi/iTA + T/), 
and similarly = 7s(,r/(so)/[r^(so) + r/fsp)]. In the 



case of an asymmetric random walk model [22| , T 



A/I 



0, 



both set of probabilities {qmi^)}m=o become Poissonian 
counting processes with 



1- 



and / 



76 



By working Eqs. ([T^ and (|20l) in a Laplace domain, 
/(^) — /(i)e~^*, it is possible to write Z+(s,i) as a 
superposition of two exponential functions scaled by the 
roots [Q{0 — 0] of the characteristic polynomial 



(21) 



From its definition, the smaller solution [22|, after the 
change of variable s )^/SI, provides 6+(A). We get 



e+(A) 



e' + Fa + F/ 



1 



1/2 



(22) 

where 9'^ = Ox/si = 7(1 - e^^/^^). This resuh, joint 
with Eq. (IT^ and the transformation ([TT|) completely 
characterize the asymptotic statistics and LDF of 

In Fig. 2(a) we plot the LDF 9 (A) for two different 
set of parameter values {sq, 7) T^, F/}. In both cases the 
symmetry ([T2|) is satisfied. The values of a were deter- 
mined from Notice that 8(A) cannot be well ap- 
proximated by a quadratic polynomial, which in turn im- 
plies the presence of strong non-Gaussian features. Figs. 
2(b) and 2(c) confirm this fact. Both LDF finr) were 
obtained numerically through the Legendre transforma- 
tion (ITlT) . The same results follows from the asymp- 
totic behavior of the exact distribution P^fJ-r) obtained 
by Fourier inversion of its generating function. Consis- 
tently, the LDF and the asymptotic exact distribution 
satisfy the FR ^TU\\ . In the insets, we plot the time de- 
pendence of the skewness 5* = (Sn^) ja^ and kurtosis 
K = {SfJ-r} /f^ ~ 3, obtained from the generating func- 
tion of P{iir), where Sfir = fJ-r — (Mt) and — (Sis^) . 
These objects also confirm the strong non-Gaussian na- 
ture of P{^J.t)- 

The functions ^(fJ-r) not only depart from a quadratic 
polynomial, but also develop a "kink," that is, and 
abrupt change around the origin. This characteris- 
tic arises when &{X) presents a plateau regime cen- 
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FIG. 2: (color online) (a) LDF e(A) [Eq. ^] for two differ- 
ent values of a. (b)-(c) LDF <p{fj,r) [Eq. ^] for a = 43. 7s"^ 
and a = 65.5s~^ respectively, (d)-(e) Realizations of fJ^rit) 
[a = 43.7s~^] for r — 0.025s and r — 0.1s respectively. In 
all cases, the dotted lines follows from approximation (|23p . 
The insets show the skewness (S) and kurtosis (K) as func- 
tion of T indicating the strongly non-Gaussian nature of /ir- 
The parameters [Eq. (|20[) ] for the curves with a = 43.7s~^ 
are so = 1.56, 7 = AOs'^ Ta = 0.54s-\ and F/ = 1.25s'\ 
For Q = 56.3s~^ are the same except Tj = 5s~^. 



tered around a/2. For the chosen parameters values, the 
curves and these special features are similar [2^ to those 
found in the numerical and experimental results of Refs. 
[2^ [23} . In the present model, a physical effect can be 
associated to these properties: the kink and the plateau 
regime are closely related with the development of inter- 
mittence in the stochastic realizations of /ir • This is the 
third main result of this contribution. 

The stochastic realizations of n'^ti'^) can be obtained 
from Eq. ([^0]) by using standard monte Carlo meth- 
ods. Alternatively, a more simple algorithm follows by 
mapping Eq. (PO)) with a renewal process |26|]. On the 
other hand, the realizations of n~^(t) can be obtained 
from a conditional scheme defined in Ref. {26| , where re- 
alizations of n^t(t) with m-events are selected as one of 
n^ti^) with probability e^""™ and discarded with proba- 
bihty (1 — 6^"°™). Statistically independent realizations 
of nf^{t) allow to generate the stochastic trajectories of 
fir [see Eqs. (g]) and (fT6)) ]. 

In Figs. 2(d) and (e) we plot two realizations of /ii-, 
corresponding to a = 43.7s~^, for two different values 
of T. Consistently with the definition ^ they fluctuates 
around — 1. For the chosen parameter values, the 
trajectories switch amongst periods of time where fir is 
positive, negative or null. As expected, for increasing r 
the fluctuations are diminished. Furthermore, for higher 
averaging times (not shown) the realizations only assume 
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FIG. 3: (color online) (a) LDF e(A) [Eq. ^] with a = 
12.5s~\ (b) LDF (fiifir) [Eq. ^]. (c) Probability distribu- 
tion of fir for r = 0.22s and r — 0.70s. (d)-(e) Realizations 
of i-irit) for r = 0.1s and r — 0.3s respectively. Inset show 
the skewness (S) and kurtosis (K) as function of r indicating 
tendency to become more Gaussian at higher time r. In all 
curves the parameters [Eq. (|20|) ] are so = 0.83, 7 = 25s~^, 
Ta = 0.3s~^, Ti = 20s^^. The dotted lines follows from ap- 
proximation (|23p . 

positive values. From Eq. (j20l) . it is immediate to deduce 
that when 7 ^ {Ta, Tj), the realizations of develop 
intermittence. This is not a surprising result. What is 
novel and not trivial, is that one can build up the com- 
plementary process n~f{t) such that their (normalized) 
subtraction [Eq. ([5])] satisfy the FR. The closeness be- 
tween the kinks of Figs. 2(b) and (c) and the intermit- 
tence shown in Figs. 2(d) and (e) is also supported by 
the results shown in Fig. 3. Here, F/ ^ F^. Therefore, 
the inactive regime is statistically (dynamically) inhib- 
ited, which in turn implies that neither n^t(t) or n^t{t) 
develops any intermittence. As shown by the plots, the 
LDF converge to quadratic polynomial functions, while 
the probability distribution becomes Gaussian. Consis- 
tently, the trajectories. Figs. 3(d) and (e), for any value 
of T do not develop the phenomenon of intermittence. 

Eqs. (Uni) and (1^^ provide an analytical expression 
for the LDF 0(A) that is very complicated. Further- 
more, the transformation Eq. pT|) . which delivers the 
LDF Lp{pLr), is only manageable numerically. Neverthe- 
less, when Ta ^ (TiT/), from Eq. ([2T|) it is simple to 
deduce the expression 

lim e+(A) = min[6l^,r/]. (23) 

This result, joint with Eqs. (fTTj) and ([TS]) allow us to ob- 
tain simple analytical expressions for the LDF 0(A) and 
(fiiHr) (see Appendix). They correspond to the dotted 
lines of Figs. 2 and 3. Evidently, they provide a very 



good fitting to the plotted curves. Nevertheless, what is 
more relevant is the possibility of getting an analytical 
description of the relation between intermittence and the 
non-Gaussian features of both LDF. 

In the limit F^ — ^ 0, both 0(A) and fdj-r) present 
points where they are not derivable functions with re- 
spect to their arguments. In particular, when intermit- 
tence develops, Tj <^ 7, the LDF (fiifir) has a hnear be- 
havior around the origin with different slopes for /j,t- ^ 
[see dotted lines (/^t ^ 0) in Figs. 2(b) an (c)]. The 
addition of the slopes of (/'(Mt) around the origin is equal 
to —a [see Eq. (IA.12P ]. The experimental data of Ref. 
[23| seem to be consistent with this relation. On the 
other hand, in the intermittence regime 0(A) develops a 
plateau regime [see dotted lines in Figs. 2(a)] with value 
T] [Eq. (jA.lip j. which is similar to that found in the 
numerical results of Ref. [1^ . 

Finally, let us remark that the discontinuities of the 
derivative of 0(A) can be read as phase transitions be- 
tween different dynamical regimes of [IE HEl • These 
features, as well as our main proposal, Eq. ([5|), suggest 
an interesting bridge between FR and equilibrium ther- 
modynamics. 



IV. CONCLUSIONS 

In conclusion, we have demonstrated that FR are al- 
ways satisfied by the subtraction of two independent 
stochastic variables whose probability distributions are 
related by a thermodynamic-likc change of measure. This 
relation leaves completely arbitrary one of the distribu- 
tions, being this degree of freedom the basis of the present 
approach. By choosing one of then as a random mod- 
ulated Poisson process, the distribution of interest de- 
velops strong non-Gaussian features, which in turn arise 
in the parameter regime where intermittence develops. 
Hence, variables such as entropy production rates may 
satisfy FR and develop intermittence. 

Depending of the parameters values of the model the 
LDF develop a rich variety of functional dependences, 
which in turn are similar to those found in recent numer- 
ical and experimental results [29| . In particular, in the 
intermittence regime the LDF of the probability distribu- 
tion and generating function are characterized by a kink 
at the origin and a plateau regime respectively. While our 
analysis does not relies on any specific (non-linear many 
body) dynamical model, it strongly suggest that inter- 
mittence may be a central ingredient in nonequilibrium 
steady states characterized by non-Gaussian fluctuations. 

Due to the wide class of systems where intermit- 
tence develops, added to the central role of FR in non- 
equilibrium states, the characterization of dynamical con- 
ditions that guaranty the coexistence of these two prop- 
erties becomes a very interesting open issue, for which 
this paper intends to contribute. On other hand, the in- 
clusion of statistical correlations that preserve the FR in 
the two-variables model as well as the analysis of these 
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ideas in the context of deterministic thermostated sys- 
tems (phase space contraction conjecture) are addi- 
tional open issues that may deserve extra analysis. 
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Appendix: LDF in the limit ^> 

Here we obtain the LDF associated to Eq. ((23|) . 6(A) 
follows from Eq. ([T3|) . Hence, ifiifir) is obtained through 
the Legendre transformation (jlip . In order to simplify 
the expressions, we write 



e(A) = e(s = x/si), 

and similarly ifinr) as 



(A.l) 



(A.2) 



where SI so{I+ - /") [Eq. The functions &{s), 

and (p{k) are the LDF of the stochastic variable kst{t) — 
n^tit) — njf.{t)]/t. Furthermore, we define the parameters 



^ -Infl-ii 
= so + lnfl- — 



and respectively 



fc+ = 76 = 7(1 - 



k = 76 * = 7- 



£1 

7 



(A.3) 
(A.4) 

(A.5) 
(A.6) 



The LDF of a Poisson process with rate 7 is denoted as 



0(s)^7(l-e-^), 
and its (symmetrized) Legendre transform as 

Finally, we introduce the function 



(firwik) = 7(1 + e-'° - (4e-^" + — 



■In 



rl 



7 



r 



(A.7) 



(A., 



(A.9) 



Both 0(s) and (p(fc) must be defined in different param- 
eter regimes. 



i) In the parameter regime 



< — < (1 -e-"«/2), 
7 



it follows s+ < s , with 

e(s) = 

and respectively 




ipp{k) - ksQ, 
P/ — ks~ , 

ri-ks+, 

fpik), 



s < s+, 
s+ < s < 
< s, 



k < ~k+, 

-k+ < fc < 0, 
< fc < 
k+ < k. 



(A.IO) 



(A.ll) 



(A.12) 



This is one of the more interesting parameter regimes, 
which in fact corresponds to the intermittence one (Fig. 
2). Notice that in the plateau regime of 6(A) {s~^6I < 
A < SIs~) it assumes the value P/. On the other hand, 
around the origin V'(Mt) has a linear behavior which dif- 
ferent slopes for fir ^ 0. This property gives rise to the 
characteristic kink shown in Figs. 2(b) and (c). From 
Eqs. (|A.2[) and (IA.12p . we deduce that the addition of 
the slopes is equal to —a [Eq. ([T^]. Furthermore, from 
Eq. (jA.12[) . it is possible to demonstrate that the ori- 
gin is the unique point at which the derivative of fifir) 
is a discontinuous function. Remarkably, in the follow- 
ing parameter regimes ^pip-r) has a continuous derivative. 
Hence, a kink related to a discontinuous derivative only 
arises in the present case, 
ii) In the parameter regime 



(1 _e-'^o/2) < £i < (1 -e-'^"), 
7 

it follows s~ < s+, with 

0is), 

Q(s)^ {0{s)+e{~s + so)~Ti, s- 
9{-s + so), 



(A.13) 



s < s , 
< s < 



and by defining the parameters 

k* = k--k+, A = P/-6'(so), 

(k* > 0), we write 



(A.14) 
(A.15) 



ipp{k) - fcso, 
(7 — fc^) — fcs+, 
(pik) = ^ ^rw{k) - A, 

(7 - k^) - k{so - s^ 
fpik), 



k < ~k^, 
-fc" < fc < -fc*, 
-k* <k < fc*, 
k* < k < kr, 
fc" < fc. 

(A.16) 

In this intermediate regime, the non-Gaussian properties 
are gradually lost. The plateau regime of 6(A) becomes 
bend. On the other hand, we remark that even when 
fifj-r) is defined by parts, its first derivative is a contin- 
uous function. 



7 



iii) In the parameter regime 



(l-e-'*") < — < 1, 
7 



(A.17) 



with s < s"*", we get 



e{s)+A, s<s-, 
Q{s) = { e{s) + e{-s + so) - 0{so), s- <s< s+, 



(-s + so) + A, 



(A.18) 



and 



(Pp{k) — ksQ + A, k < —k , 

(7 - fc") - A;s+ + A, -k-<k<-k*, 
frw{k), —k*<k<k*, 
(7 - fc-) - fc(so - s+) + A, k*<k<k-, 
(Pp{k) + A, fc" < k. 

(A.19) 

Fig. 3 falls in this regime, where both LDF approach 
quadratic functions. In fact, while the derivative of 6(A) 



presents smooth discontinuities they are far beyond of 
the origin [lim(rj/7)^i = ±00]. As in the previous 
case, (fii^ir) has a continuous derivative. Its dependence is 
mainly defined by the function iprw{k) [lim(r//7)-!-i k* — 
00]. 

iv) In the parameter regime 



1 < < CXD, 

7 



for any value of s, we get 



(A.20) 



(A.21) 



Therefore, in this approximation 0(A) correspond to the 
LDF of an asymmetric random walk. Furthermore, it 
follows <f>{k) = iprw{k) [Eq. (|A.9p ]. As in the previ- 
ous regime, the rate fluctuations do not induce any non- 
Gaussian feature because their characteristic time is com- 
parable to that of the individual events. 
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